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Abstract 


We examine the behavior of correlation functions for a massless scalar field in de Sitter space with a quartic interaction. 
We find that two-loop corrections are relevant, and the resummation of these corrections generates a complicated structure 
whereby high momentum modes stay massless while low momentum modes develop a dynamical mass. 
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I. INTRODUCTION 


The goal of this paper is to address the following question: what do the correlation functions look like for a massless 
field in de Sitter space with a </> 4 interaction? 

This question is of more than merely formal interest; it has direct connections to issues in inflationary theory. 
The inflationary geometry has many similarities to de Sitter space (especially in the earlier stages of inflation). 
Furthermore, scalar fields play a fundamental role in inflation fli], and correlation functions of these fields are recorded 
in the cosmic microwave background. 

Current observations of the CMB from the WMAP experiment Q are currently almost exactly consistent with a 
scale invariant power spectrum. However, it would be surprising if the inflaton had no interactions whatsoever, and 
such interactions are expected to produce deviations from exact scale invariance Q. Several experiments have been 
searching for such nongaussianities. For comparison of these observations to theory, it is important to be able to 
calculate correlation functions of these scalar fields in the inflationary background. 

Unfortunately, the scalar fields in inflation are typically light (to -C H), and for these small masses, it turns out 
that the loop corrections are important. In fact, loop corrections involving massless fields are infrared divergent jdj-Q] ■ 
To compare theoretical predictions to observations, it is crucial to be able to treat these divergences rigorously. This 
has led to a lot of work using many different approaches to try and resum these divergences formally both for scalars 
and gravitons [z-44|. In some approaches the fact that one loop corrections are divergent for a massless field results 
in a self-consistent mass being developed. In other approaches, the propagator becomes time dependent, breaking de 
Sitter invariance, and in yet other approaches particle production leads again to an induced mass for the field. All 
these techniques qualitatively lead to the result that the massless fields pick up a nonperturbative mass of order y/\H 
once interactions are included. 

However, one would expect that a short wavelength mode of the scalar field should be completely unaffected by the 
de Sitter curvature (in much the same way as current lab experiments do not take into account the curvature of the 
universe). A long wavelength mode, on the other hand, can plausibly be affected by the curvature and may pick up a 
dynamical mass (indeed, this must happen in order to prevent divergences from occurring.). In short, we expect any 
corrections to be momentum dependent, with large effects only for long wavelength modes. However, the loop effects 
found in the various approaches above often seem to affect all wavelengths equally, and usually the short wavelength 
modes acquire a mass. 

The reason for this is that the analyses in most of the above papers typically focus on one-loop corrections. 
However, in c/> 4 theory, the one-loop corrections have a very non-generic feature; they are independent of the external 
momenta. This leads to the surprising effect detailed in the previous paragraph. This feature suggests that the one- 
loop corrections may not encode the relevant physics; indeed, all one-loop corrections in (jfi theory can be completely 
canceled by the addition of a mass counterterm. 

This further suggests that the resolution of the divergences involves a study of the two-loop corrections in this 
theory. Unlike the one-loop corrections, the 2-loop corrections are not independent of the external momenta. The 
long wavelength modes will then receive mass corrections which are different from the corrections to the short distance 
modes; this then allows us to cancel off the mass for the short distance modes using the counterterm, while maintaining 
a nonzero mode for the long distance modes which can cancel the IR divergence. 

Our focus in this paper is to analyze in detail how these 2-loop corrections affect the correlation functions. Unsur¬ 
prisingly, they are much more computationally difficult than the one-loop corrections, and so we will be unable to 
find closed form analytic results. However, many qualitative features can be extracted. In particular, we show that 
one loop corrections are indeed canceled by counterterms, while two-loop corrections produce a momentum dependent 
mass which is zero for the short distance modes. Hence the long wavelength modes acquire a dynamical mass, while 
the short wavelength modes remain massless, as expected. 

To see this effect, which is nonperturbative, we will need to resum the two loop corrections. Technically, we resum 
a subset of the two-loop diagrams which are expected to produce the leading infrared divergences. This resummation 
will be found to produce a momentum dependent effective mass. The scaling of the masses turns out to be very 
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different from previous calculations in the literature. 

As we shall discuss below, the difference between our results and the results in the literature appears to hinge 
on a different definition of the term ’massless scalar’ (most importantly, reference 4l| employs techniques extremely 
similar to this paper, but with a different definition of the massless limit). To define what we mean by this term, we 
use the general feature that the de Sitter propagator, on scales smaller than the Hubble scale, should locally resemble 
a Minkowski space propagator. In our definition, a massless field is one which for shorter length scales (i.e. distances 
parametrically smaller than the Hubble scale) has a propagator similar to a massless Minkowski space propagator. 
This is not the only possible choice (reflecting the difference with other definitions in the literature), but is a well 
defined one. The more naive definition would have been that the propagator should approach the propagator for a 
massless field at low momenta, but this appears to run into divergences. 


II. REVIEW OF THE IN-IN FORMALISM 


We begin by reviewing the in-in formalism for ^ 4 
already been derived and presented elsewhere (e.g. 
further details. 

We will take the metric of de Sitter space to be 


theory. This section deals with the free field. These rules have 
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30[); we refer the interested reader to these papers for 


ds 2 


H 2 i 


dr 2 — ^2 d x i 


(1) 


We will also consider a scalar field propagating in this geometry; the field will have mass m and a quartic interaction. 
The Lagrangian is then 


m = v~a 


^ v d^(t)d v (t> - ^ m 2 (j) 2 



( 2 ) 


We note that we will be focusing on very light fields. To leading order, therefore, we will set m 2 = 0. Furthermore, 
we will have to work in the in-in formalism. In this formalism, the number of fields is doubled; for the scalar field, we 
have the fields (j >+, . We additionally define 


&C = 2 ^+ + &-) 0A =<!>+ ~<t>- 


( 3 ) 


The in-in Lagrangian is then defined as £ = C(<f>+) — £(^_). 

There are four propagators corresponding to the two sets of scalar fields. The Keldysh propagator is denoted F 
and is defined as 


iF(x,y) = {(f)c{x)(t>c{.y)) 

In addition, we have the advanced and retarded propagators 

G R {x,y) = {<t>c{x)<l>A(y)) = iO{x 0 - y°)((<t>(x)<t)(y)) - (<?%)#*0» 
G A {x,y) = G R (y, x) 


( 4 ) 


( 5 ) 

( 6 ) 


The fourth propagator (4>A(x)4>A(y)) is identically zero. 

In Feynman diagrams, (j>c is denoted by a solid line, and 4>a by a dashed line. The propagators (after a Fourier 
transform) are denoted by 


= F(k, Ti, 72) 
= - iG R (k,T 1 ,T 2 ) 


( 7 ) 

( 8 ) 


There are two vertices in <^> 4 theory. 
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\ / 


( 9 ) 




( = — iXa 4 {r) 


( 10 ) 

We can calculate these propagators in the free field limit i.e. when A = 0. For free fields, the propagators can be 
found to be (we will denote the free field limit by a 0 superscript) jill 


■kH 


-F—2 (k, Ti, 72 ) = —— (rir2) 3/2 J Re(iF 1/ (-/cTi)iJ*(-fcr2)) 


G^ ( 2 0) (fc, Ti,T 2 ) = —9{t\ - T 2 )^—(TlT2) 3/2 Im(H l ,(-kT 1 )H*(-kT2)) 


7 tH 


( 11 ) 

( 12 ) 


where v 2 = \ — jj?- 

There are several limits of the free field F-propagator that are of interest. The first is the zero mass limit. For zero 
mass, we have 


F^ = 0 (k, T i, r 2) = 7773 [(1 + fc 2 nr 2 ) cos(fc(n - t 2 )) + fc(n - r 2 ) sin(fc(n - r 2 ))] 


H 


2 k 3 


G^ 2 > l 0 (k,T 1 ,T 2 ) = 0(n - t 2 )1^-[(1 + fc 2 Tir 2 ) sin(fc(ri - r 2 )) - k{r 1 - r 2 )cos(fc(ri - r 2 ))] 

1 

W 


H 2 


(13) 

(14) 


^ m 2 = 0 V'"? ' -L’ ' 3 

For small k the Keldysh propagator goes as A 

P (o) n \ # 2 

^ m a=oA,n,T 2 ) - 

This leads to divergences in loop diagrams. For finite masses, the low momentum behavior of the propagator is 
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expressible in terms of = e<l and is found to be 

H 2 


(15) 


F^{k,T!,T 2 ) ~ ^{k 2 T lT2 Y 


(16) 


which regulates the infrared divergences found at m 2 = 0. 


III. ONE-LOOP CORRECTIONS 

We first briefly consider the one-loop corrections to the massless theory. The relevant diagrams are shown below. 
The third and fourth diagrams represent counterterms. 



♦ 


♦ 


FIG. 1: One loop corrections 


In a perturbation expansion, the first diagram is of the form 



(t-3)G (0) (/c,ti,t 3 ) 


d 3 qF^(q,T 3 : r 3 ) 


F^(k,T 3 ,T 2 ) 


(17) 


For to 2 = 0, F( 0 )(g, t 3 ,t 3 ) ~ 4j- at small q, and so the integral is divergent. The massless theory therefore does 
not seem to have a good perturbation expansion. However, the infinite correction can be canceled by a suitable 
choice of the counterterms. Since the loop diagram is independent of external momenta, canceling the mass at high 
energies also cancels the mass at zero momentum. This will again result in a massless zero mode which will lead to 
IR divergences. We must therefore go beyond one loop to see whether the divergences can be resolved. 
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IV. TWO-LOOP CORRECTIONS 


The two loop correction has the generic form of a sunrise diagram. There are in fact three separate two-loop 
corrections which occur in this theory, which we denote £ 1 , 2 , 3 - 

The first set of corrections is shown in Figure 2. These corrections can be written explicitly as 




FIG. 2: Corrections contributing to Ei 


Si(fc,r 1 ,r 2 ) = *A 2 a 4 (r 1 )a 4 (r 2 ) J d 3 k ! J d 3 k 2 F(ki,r 1 ,T 2 )F(k 2 ,T 1 ,T 2 )G R (k - h - fc 2 ,r 1 ,r 2 ) 
-?Y a4 ( T i) a4 ( T 2) J J d 3 k 2 G R (ki,Ti,T 2 )G R (k 2 ,T 1 ,T 2 )G R {k-ki-k 2 l Ti 1 T 2 ) 


(18) 


This loop can be thought of as inducing a new vertex corresponding to a (nonlocal) term in the Lagrangian 
proportional to ■ 

There is an analogous vertex where we reverse the time direction; this diagram is the mirror image of the diagram 
above. This mirrored vertex will be denoted £ 2 and can be thought of as inducing a new vertex corresponding to a 
(nonlocal) term in the Lagrangian proportional to 

Finally, there is a third set of loop corrections which induce a new vertex corresponding to a (nonlocal) term in the 
Lagrangian proportional to This set of diagrams will be denoted £ 3 , and can be written as 

£ 3 (fc,ri,r 2 ) = ^-a 4 (n)a 4 (r 2 ) J d 3 ki J d 3 k 2 G R (k 1 ,T 1 ,T 2 )G R (k 2 ,T 1 ,T 2 )F(k - ki - fc 2 ,n,r 2 ) (19) 

+ ja 4 (n)a 4 (r 2 ) J d 3 k 1 J d 3 k 2 G A {k 1 ,T 1 ,T 2 )G A (k 2 ,T 1 ,T 2 )F(k - ki - fc 2 ,n,T 2 ) (20) 

-A 2 a 4 (ri)a 4 (r 2 ) J d 3 k x J d 3 k 2 F(ki,T 1 ,T 2 )F(k 2 ,T 1 ,T 2 )F(k - h - k 2 ,T\,T 2 ) (21) 


For compactness, we shall replace all loops £; in Feynman diagrams by a insertion labeled with the appropriate 
index i. Therefore for example, we shall indicate a £3 insertion by a dot with a 3 subscript as shown in figure 3. 



FIG. 3: A condensed notation for the loop integrals 


V. THE HIGH MOMENTUM REGION 

These diagrams are difficult to compute exactly. We are, however, interested in the most divergent pieces of these 
diagrams; this will allow us to make several approximations to obtain an estimate of these loop corrections. We will 
focus on £1 at the beginning, and apply the same approximations to the other integrals afterward. 
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We first note that the IR divergences primarily come from regions of the loop integral where one of the F propagators 
has a small momentum flowing through it. Indeed, if we were to replace F by F ^ in the above integrals, they would 
all be infrared divergent because F^°\k, t\, 72 ) ~ -p-. This IR divergence is presumably regulated when the full F 
propagator is used. We therefore conclude that the low momentum behavior of the full F propagator is significantly 
different from the tree level propagator F°. This modification by the higher order corrections will moderate the 
infrared behavior and make the integrals finite. Now, although the F propagator is modified in such a way that the 
integrals are finite, the divergences must reappear if we take the interaction A —> 0. It must therefore be the case 
that we get a large contribution (i.e. parametrically larger as a function of A) from the region of integration where 
the momenta flowing through an F propagator is small. This further indicates that we can keep only the terms with 
the largest number of F propagators, which are expected to have the largest contributions. 

On the other hand, the G propagator does not have infrared divergences when the momentum flowing through it is 
small. This means that there is no requirement for the corrections to modify it significantly at low momenta. We will 
assume the corrections to G are small (i.e. suppressed by powers of A) which will be confirmed by later calculations. 
We will therefore (to leading order) replace G by G° in the above expressions. 

Furthermore, the integration over momenta is expected to be dominated by the regions when both F propagators 
have small momenta in their arguments. In particular, we would expect the dominant part of the integral for £1 to 
come from the integration region where both internal momenta fci, fc 2 are much smaller than the external momentum 
k i.e. fci, fc 2 < k. The integral then simplifies to 

£i(fc,n,T 2 ) ~ ia 4 (T 1 )a 4 {T 2 )G R ^ ) (k,T 1 ,T 2 ) x [c(fc, n, r 2 )] 2 (22) 


where we have defined 

fk 

c(k, ti, 72) = A / d 3 k 1 F{k 1 ,Ti,T 2 ) 
Jo 


(23) 


The next point to note is that the leading term in the c(fc, ti, t 2 ) integral above is expected to be time independent. 
The integral diverges if F is replaced by the tree level F° propagator, which for small k leads to an integral of the 
form 4^r■ The divergent term is time independent. When the divergence is regulated, the integral will be large but 
finite; it is reasonable to expect that the leading value (in an expansion in A) will continue to be time independent. 
We will defined this leading term to be m\{k). 

With these approximations, we are able to simplify the £1 integral to 


£i(fc,r 1 ,T 2 ) ~ 7:a 4 (ri)a 4 (r 2 )?ni(fc)G K(0) (/c,ri,T 2 ) 


(24) 


We can now resum the two loop corrections. For the retarded propagator, the corrections that we wish to sum are 
of the form shown in Figure 4. 


— ■•f—.? - — 

FIG. 4: Corrections to the retarded propagator 


These lead to an equation for the full propagator 

G R (k,T 1 ,T 2 ) = G R ^°\k,T 1 ,T 2 ) (25) 

+m\{k) J (jjV^)4 GR{0) w r 3)G fl(0) (fc, 73 , r 4 )G fl(0) (fc, r 4 , r 2 ) (26) 

+m?(t) / (gr^>- (Hr))- (g^)< (27) 

xG fl(0) (fc, 75 , T 6 )G fl(0) (fc, T 6 , T 2 ) (28) 

+ ... (29) 
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We can perform the sum by comparing it to the expansion when a mass term is treated as a perturbation. A 
massive field has the exact retarded propagator G m 2 (k, t\, r 2 ). On the other hand, in the same theory, the mass 
term —m 2 (f>c(f >a can be treated as a perturbation to the massless theory. We can match the exact solution to the 
perturbation expansion; this yields 


Gy! ( 2 0) 0, n, r 2 ) = G^ ( 2 °2 0 (fc, ti , r 2 ) 


— m 


dr 3 ^r(o) 


Gyi 0 {k, n , T 3 )G^l 0 (k, r 3 , r 2 ) 


(-Ht 3 ) 4 m2=0V 

4 /' dr 3 dr 4 i?(o) ,, 0) ^.R(o) /, \ 

+ J (ffr 3 ) 4 (gr 4 ) 4 o( fc >n,r 3 )G m V 2 =0 (fc, t 3 ,T 4 )G m 2 _ 0 (fc,r 4 ,r 2 ) 
+ 


(30) 

(31) 

(32) 

(33) 


Comparing this to the previous expansion involving Si, we find 


G fl (fc,Ti,T 2 ) 


2 l 


(fc,Ti,T 2 ) + G 


R( 0 ) 

— (fc) 


(k,T 1 ,T 2 )) 


(34) 


-■j— ? - »—t - r—t— 


FIG. 5: Corrections to the Keldysh propagator 

Similarly, the F propagator receives contributions of the form shown in Figure 5. Just like the the retarded 
propagator, we find that the contribution of these terms is ^(F m 2 ^(k,Ti,T 2 ) + F_ m 2 ^(k, ly, r 2 )). 

We also need to consider contributions which involve a S 3 vertex insertion. These are all of the form shown in Fig. 
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-•? • -3- ••• f-T" 


FIG. 6: Other corrections to the Keldysh propagator 


Comparing to the expansion for the retarded propagator, we find that these contributions can be written as 

j -^^Aj^j4 GR (k,T 1 ,T 3 )^ 3 {k,T 3 ,T4)G A (k,T 4 ,T2) (35) 

We shall treat these diagrams in perturbation theory; that is, we shall not include them in the resummation. We will 
hence ignore these terms for the remainder of the discussion. It would be interesting to examine the effects of these 
diagrams; we shall leave this for future work. 

Summarizing our results so far, we have found that the resummed Keldysh propagator is of the form 

F(k,T U T 2 ) = i(K m 2 (fe) (fc,n,r 2 )+ F_ m 2 (fe) (fc,n,r 2 )) (36) 

where m\(k) is the leading term in the integral (1231) . Note that this is not the propagator for a massive field. 


VI. THE LOW MOMENTUM REGION 

Unfortunately, our result above cannot hold for all external momenta. For if it did, (1361) would be the exact 
solution for the propagator. But now the term in the propagator F_ m 2^ (k, r 4 , r 2 ) grows in the infrared even faster 
than Fo(fc, r 4 , r 2 ). From E3l) . we would find m\{k ) to be divergent. This implies that the solution is not self consistent. 
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One possibility that there is simply no solution; that there is a deep inconsistency that does not allow a propagator 
to exist when the one-loop corrections are canceled by the counterterms. While this possibility cannot be neglected, it 
would be surprising, for the reasons outlined in the introduction. We therefore conclude that some other corrections 
must become relevant at low momentum. The consequence is that there must be a crossover momentum ko below 
which our solution (13611 breaks down. 

We can trace this breakdown back to the original simplifications we made for the integral. We found that the 
leading contribution to £1 was 

Si(fc,ri,r 2 ) ~ iA 2 a 4 (ri)a 4 (r 2 ) J d 3 ki J d 3 k 2 F(ki, n, r 2 )F(k 2 , n, t 2 )G r(0) {k - k x - k 2 , n, r 2 ) (37) 

We assumed that the main contribution came from the region of integration where the momenta k\, k 2 were very 
small, and in particular much smaller than the external momentum k. Under this assumption, the integral was found 
to be proportional to [c{k, ti, t 2 )] 2 where c(fc,ri,r 2 ) was the integral of F(k\,Ti,T 2 ) up to momentum k. But for 
small external momenta k, the last factor necessarily goes to zero. It therefore cannot be dominant for small external 
momenta. Once again we are led to the existence of a crossover momentum ko below which the solution (13611 is not 
valid. 

On the other hand, for larger momenta (k above ko), our resummation in the previous section is still valid. We 
then expect the form of the propagator at these momenta to be of the form (l36l) . 

For external momenta which are smaller than ko, we can approximate the integral for £i by its value at k = 0: 

/•OO 

£i(A;,ti,t 2 ) ~ iA 2 a 4 (ri)a 4 (r 2 ) / d 3 k 1 d 3 k 2 F(ki, n, r 2 )F(k 2 , n, T 2 )G R {-k\ - k 2 , n, r 2 ) (38) 

J fci,/c2=0 

This is independent of k, and hence produces a correction that is dominant for low momenta. These corrections will 
therefore produce an effect similar to a dynamical mass, and in particular, they will cut off the infrared divergences. 
The time dependence is more complicated than a simple mass insertion, and so we therefore make the ansatz 


F{k,Ti,T 2 ) 


^{F m l(k)(k,Ti,T 2 ) +F_ m 2 (fe) (fc,n,r 2 )) k > k 0 
f(Ti,T 2 )^{k 2 Y k<k 0 


(39) 


where dgy = e, Yfp = e> j an d we have parametrized the time dependence at low momenta by the function /. 

We now can attempt to find the scaling of the various parameters of our ansatz with A. Comparing the equations 
(l23l) . which contains one integral over A F, and (l38l) . which contains two integrals over A F, we expect m 2 ~ ( m f) 2 , 
and thus that m\ is parametrically smaller than m 2 . In turn, this means that the integrals in these equations are 
dominated by the low momenta regions k < ko- We then find 

ml(k) oc A ( fc o)“ £ (^ 7 ) m l ~ A 2 (fc^)" 2e (-^) (40) 

Since e' is parametrically small in A, we can set to leading order (fc^) -26 = 1. We then find m 2 ~ A 2 / 3 and m\ ~ A 1 / 3 . 

Our conclusion for the low momentum modes is then that their propagation is modified in a manner similar to a 
dynamical mass. This dynamical mass has a scaling proportional to A 1 / 3 and cuts off the infrared divergences. 


VII. DISCUSSION AND CONCLUSION 

We have discussed how a quartic interaction would modify the propagator of a massless field in de Sitter space. We 
have found a complicated result; there is a crossover momentum scale with very different behavior of the propagator 
above and below the crossover. For long wavelengths, the modes behave as if they develop a dynamical mass. For 
shorter wavelengths, surprisingly, the propagator is the sum of a massive propagator and a propagator for a tacliyonic 
field. Our main result is encapsulated in equation (1391) . which is our result for the Keldysh propagator once interactions 
are included. 


The apparently surprising behavior at shorter momenta has a natural explanation from the masslessness of the 
scalar. Consider the behavior of the free Keldysh propagator for iq = r 2 = r with kr large. The free propagator has 
the expansion 


F ^2 {k, ri, r 2 ) = ^7 — [1 + 


m 


>) + •••] 


(41) 


2k ' k 2 r 2 2 k 2 H 2 r 2 ' 

We can therefore define the mass of the particle by looking at the deformation away from the massless propagator i.e. 


m 2 = lim 4 k 3 [F(k, n, r 2 ) - F^L 0 (fc, n, r 2 )] 


(42) 


For the expression (136[) . we find that the mass is zero, thereby justifying calling this a massless scalar. At this order 
in the loop expansion, the high momentum modes in our calculation do not obtain a mass correction. Note that the 
tachyonic piece was necessary for this to happen. 

The behavior for longer wavelengths is more in line with expectations. We have found that momenta with longer 
wavelengths behave as if they have a mass, which is similar to results in the literature. Quantitatively, however, our 
results are different; we have found a dynamical mass which scales as A 1//3 rather than the usual scaling of A 1 11 / 4 found 
by previous authors. This is because the mass in our case is generated at two loops rather than one-loop. 

Further corrections to the propagator and other correlation functions can be calculated in perturbation theory. 
Infrared divergences are now absent and so the perturbation expansion should make sense. We note that any loop 
involving an F propagator is enhanced by a factor ~ A _2//3 . This indicates that the perturbation expansion is in 
powers of A 1 / 3 rather than A. 

There are several open issues that still need to be addressed. In particular, we have argued for a crossover momen¬ 
tum, but we have not found its magnitude (even as a scaling in A). A calculation to next order in perturbation theory 
might shed light on this issue. We shall leave this and other questions for future work. 
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